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We analyze the spectral properties of the quark propagator above the critical temperature for the 
deconfinement phase transition in quenched lattice QCD using clover improved Wilson fermions. 
The bare quark mass dependence of the quark spectral function is analyzed by varying the hopping 
parameter K in Landau gauge. We assume a two-pole structure for the quark spectral function, 
which is numerically found to work quite well for any value of K. It is shown that in the chiral 
limit the quark spectral function has two collective modes that correspond to the normal and 
plasmino excitations, while it is dominated by a single-pole structure when the bare quark mass 
becomes large. 
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1. Introduction 

Hadronic matter undergoes a phase transition to a deconfined phase at nonzero temperature. 
The properties of the phase above the critical temperature, Tc, acquires much interest both from the 
experimental and theoretical point of view. In order to understand the structure of the matter in this 
region, it is desirable to identify the basic degrees of freedom of the system and their quasi-particle 
properties. In the present study, we analyze dynamical properties of quarks above Tc in quenched 
lattice QCD in Landau gauge |[T]]. 

At asymptotically high temperatures one can calculate the quark propagator using perturbative 
techniques. It is known that the collective excitations of quarks in this limit develop a mass gap 
(thermal mass) that is proportional io gT |^ ||]. Here g and T denote the gauge coupling and 
temperature, respectively. Moreover, in this limit the number of poles in the quark propagator is 
doubled. In addition to the normal modes, which reduce to poles in the free particle propagator, 
plasmino modes appear. 

In order to understand the origin of the plasmino mode in the quark propagator in the high 
temperature limit, it is instructive to consider the quark propagator at intermediate temperature 
[0, ^]. In the temperature dependence of the spectral function for fermions with scalar mass m 
has been considered in QED and a Yukawa model. In these models the spectral function at zero 
temperature has two poles at energies ft) = itm, while in the high temperature limit, T/m — > oo, it 
approaches the propagator having four poles. The one-loop calculation performed in [Q] clearly 
showed that the two limiting forms of the spectral function are connected continuously; in addition 
to the normal quasi-particle peak a peak corresponding to the plasmino gradually appears in the 
spectral function and becomes larger with increasing temperature, [Q]. In the present study, we 
analyze the quark propagator in quenched lattice QCD at two values of the temperature, T = \.5Tc 
and 3Tc, as a function of the bare quark mass. To simplify the present analysis, all our calculations 
have been performed for zero momentum. 



2. Quark Spectral Function 

The dynamical properties of quarks at zero momentum are encoded in the quark spectral func- 
tion p(ftj) which is related to the Euclidean correlation function 

S{x) = i| d\(Py{^f{x,xmO,y)), (2.1) 

through an integral equation 

5(T) = j_Jco-^;^j^^^-^p{co), (2.2) 

with the quark field the spatial volume V , and the imaginary time T which is restricted to 
< T < 1/r. The Dirac structure of p {(o) is decomposed as 

p(ftj) = Po(ftj)/' + Ps(w) 

= p+(a))A+y' + p_(ftj)A_/', (2.3) 
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Figure 1: The spectral function p+{co) in the Yukawa model with massive fermion and massless boson for 
various values of m/T [0]. 
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Table 1: Simulation parameters 



with projection operators A± = (1 it 7°)/2. The charge conjugation symmetry leads to po(co) = 
Po(-ft)), Ps(ft>) = -Ps(-ft)), and p+(ft)) = p_(-G)) = po(ft)) + Ps(ft)) iQl- In the following anal- 
ysis, we concentrate on a determination of p±{o)) instead of po,s(ft)), because excitation prop- 
erties of quarks are more apparent in these channels. The spectral functions p±{o)) are neither 
even nor odd functions. In the chiral limit, however, ps vanishes and p±(ft)) become even func- 
tions. In analogy to Eq. we introduce the decomposition of the correlation function S{t) as 
S(t) = 5+(t)A+7" + 5_(t)A_7", where S± are related through 5+(t) = 5_(j3 - t). 

For free quarks with scalar mass m the spectral functions, p±{(o) = n5{(0^m), have quark 
and anti-quark poles at ft) = itm, respectively. In the high temperature limit, additional poles, 
corresponding to the plasmino, appear at negative energy for p+(ft)) and positive energy for P-(ft)) 
[0, ^. In Fig. [l| we show P+(g)) in the Yukawa model with a massless boson at finite T for various 
values of m/T [^, One sees that the shape of p+{(o) changes continuously between these two 
limits i&mjT varies. 



To extract the spectral function p+(G)) from 5(t) using Eq. (p^, we assume that P+(g)) can 
be described by a two-pole ansatz. 



p+(a)) =Zi5(a)-£i)+Z25(a) + £2), 



(2.4) 



where the residues Zi 2 and energies E\,^ > have to be determined from a fit to 5+(t). The poles 
at ft) = Ei,—E2 correspond to the normal and plasmino modes, respectively [ffl]. 
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Figure 2: The lattice correlation function 5'4'" (t) at T = 3Tc for the lattice of size 64^ x 16 with various 
values of k, and the fitting result with the ansatz Eq. (2.4). 



The correlation function S{t) has been calculated at two values of the temperature, T = \ .5Tc 
and STc, in quenched QCD using non-perturbatively improved clover Wilson fermions ^. To 
control the dependence of our results on the finite lattice volume, A^^ x A^^, and lattice spacing, a, 
we analyze the quark propagator on lattices of three different sizes. The gauge field ensembles 
used for this analysis have been generated and used previously by the Bielefeld group to study 
screening masses and spectral functions [Q]. The different simulation parameters are summarized 
in Table [l] For each lattice size, 51 configurations have been analyzed. Quark propagators 
have been calculated after fixing each gauge field configuration to Landau gauge. For this we used 
a conventional minimization algorithm with a stopping criterion, {\/3)tr\d^A^\^ < 10^'^ In the 
Wilson fermion formulation the bare mass, mo, is related to the hopping parameter K, through 
the standard relation niQ = (1/k"— /{2a), where Kc denotes the critical hopping parameter 
corresponding to the chiral limit, or vanishing quark mass. 



3. Numerical Results 

In Fig. ^ we show the numerical results for 5!^"' (t) for several values of K calculated on a 
lattice of size 64^ X 16 at r = 3T,. One sees that the shape of S*!"'' (t) approaches that of a single 
exponential function for smaller k:, while it becomes symmetric as K approaches Kc. In the vicinity 
of the wall source, i.e. at small and large T, we see deviations from this generic picture which 
can be attributed to distortion effects arising from the presence of the source. We thus exclude 



points with t < T„„„ and — r < T„„„ from our fits to the ansatz given in Eq. ( |2.4[ ). The resulting 
correlation functions obtained from correlated fits with Tmm = 3 are shown in Fig. ^. One sees 
that S^f-^- is well reproduced by our fitting ansatz^; the ;^^/dof of our fits is between 2 and 3 at 



'We also checked that fits based only on a single pole ansatz lead to unacceptable large ;f^/dof. 
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1.335 ^ K" < 1.34, while it gradually increases as K becomes smaller than K = 1.335. A similar 
behavior is also observed for our other lattice sizes [^]. 

In Fig. ^, we show the dependence of £"12 and Z2/(Zi +Z2) on the bare quark mass mo for 
T = 1.5rc and 37^. The results have been obtained from two-pole fits on lattices of size 64^ x 16. 
Errorbars have been estimated from a Jackknife analysis. The dotted line in this figure denotes 
the pole mass determined from the bare lattice mass. The figure shows that the ratio Z2/ (Zi +Z2) 
becomes larger with decreasing mo and eventually reaches 0.5. The hopping parameters satisfying 
Zi = Z2 are k'^ = 0.133974(10) for T = 37^ and = 0.134991(9) for T = l.5T„ which are consis- 
tent with the values for jq. given in Table |l|. The latter had been obtained in from a fit to critical 
hopping parameters determined in [||] from the vanishing of the isovector axial current. The nu- 
merical results obtained on 64^ x 16 lattices show that Ei and E2 are equal within statistical errors 
at K" = k'^. The spectral function p+{(o) thus becomes an even function at this point; the quark prop- 
agator becomes chirally symmetric despite the presence of a thermal mass, mj = Ei = E2. From 
Fig. ^, one also finds that the ratio mj /T is insensitive to T in the temperature range analyzed in 
this work, while it is slightly larger for lower T . 

As mo becomes larger, Z2/ (Zi +Z2) decreases and p+(«) is eventually dominated by a single- 
pole. One sees that E\ has a minimum at mo > 0, while E2 is an increasing function of mo. In the 
one-loop approximation, the peaks in p+(ft)) corresponding to E\ (£"2) are monotonically increas- 
ing (decreasing) function of m^/T [Q, |^. The quark mass dependence of poles found here thus 
is qualitatively different from the perturbative result. We find, however, that slope of E2 as func- 
tion of mo/r decreases with increasing T . This may suggest that the perturbative behavior could 
eventually be recovered at much larger temperatures. 

In order to check the dependence of our results on the lattice spacing and finite volume, we 
analyzed the quark propagator at T = JiTc for three different lattice sizes. Results for E\ and E2 
are shown in Fig. ^. Comparing the results obtained on lattices with different lattice cut-off, a, but 
same physical volume, i.e. 64^ x 16 and 48'' x 12, one sees that any possible cut-off dependence 
is statistically not significant in our analysis. On the other hand we find a clear dependence of the 
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Figure 4: The bare quark mass dependence of parameters £1, £2 at T = 37^ for lattices of size 64^ X 16, 
48^ X 16 and 48^ x 12. 
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Figure 5: The bare quark mass dependence of fitting parameters £12 and Z^j (Z\ +Z2) at T = 1.57^ and 
for lattice 64^ x 16. 



quark energy levels on the spatial volume; when comparing lattices with aspect ratio Na/Nx = 3 
and 4 we find that the energy levels, E\^2, drop significantly. A similar behavior is observed also at 
T = l.5Tc. 

The presence of a strong volume dependence of the quark propagator is not unexpected. In 
fact, the thermal quark mass arises as collective effect of low momentum gluons; gluons at the 
soft scale p ^ gT play a crucial role to give rise to the thermal mass at high temperatures 
However, on lattices with given aspect ratio Na/Nt low momentum gluons are cut-off. The low- 
est non-vanishing gluon momentum is, Pmin/T = 2n{Nx/Na), which still is larger than unity on 
lattices with aspect ratio Na/Nt = 4. The situation may, nonetheless, be somewhat better in the 
temperature range explored here as the temperature dependent coupling g{T) is larger than unity. 
An analysis of quark spectral functions on lattices with even larger spatial volume is needed in 
the future to properly control effects of small momenta. We attempted to estimate the thermal 
mass in the V —>-oo limit by extrapolating the results obtained for two different volumina. Defining 
mj = {Z\E\ +Z2E2) / {Z\ +Z2)|k:=k;, and assuming the volume dependence of my as mjiNx/Na) = 
mr(0)exp(A^3/^3)^ ^^^^^^ mj(0)/T = 0.771(18) for T = 37^ and mriO)/! = 0.800(15) for 
T = l.5Tc. This suggests that finite volume effects may still be of the order of 15% in our current 
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analysis of mj jT (see Fig. ||). Despite these problems, our result clearly shows that light quarks 
near but above have a mass gap that is of collective nature similar to that in the perturbative 
regime. 

4. Summary 

In this study, we analyzed the quark spectral function at zero momentum for T = 1.57^ and STc 
as functions of bare quark mass mo in quenched lattice QCD with Landau gauge fixing. We found 
that the two-pole approximation for p+(G)) well reproduces the behavior of the lattice correlation 
function. It is argued that the chiral symmetry of the quark propagator is restored at the critical 
value of K and the shape of the spectral function at this point takes a similar form as in the high 
temperature limit having normal and plasmino modes with thermal mass mj. Meanwhile, P+(g)) 
approaches a single-pole structure as mo is increased, as one can naturally deduce intuitively. The 
non-perturbative nature of thermal gauge fields is reflected in the behavior of poles as functions 
of mo, which is qualitatively different from the perturbative result [Q]. We also note that the ratio 
mj /T decreases slightly with increasing T , which is expected to happen at high temperature where 
mj /T should be proportional to a running coupling g{T). Although results on the quark propagator 
are gauge dependent, we expect that our results for its poles suffer less from gauge dependence, 
because the success of the pole approximation for p+{(o) indicates that the quark propagator has 
dynamical poles near the real axis, which are gauge independent quantities. 

In the present study, we analyzed the quark spectral function in the quenched approximation. 
Although this approximation includes the leading contribution in the high temperature limit 
and thus is valid at sufficiently high T , the validity of this approximation near Tc is nontrivial. For 
example, screening of gluons due to the polarization of the vacuum with virtual quark antiquark 
pairs is neglected in this approximation. The coupling to possible mesonic excitations, which may 
cause interesting effects in the spectral properties of the quark [||], are not incorporated, either. The 
comparison of the quark propagator between quenched and full lattice simulations would tell us 
the strength of these effects near Tc. 
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